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Abstract 

In order to test gravitation in the Solar System, it is necessary to improve the orbit restitution 
of interplanetary spacecrafts. The addition of an accelerometer on board is a major step toward this 
goal because this instrument measures the non-gravitational acceleration of the spacecraft. It must 
be able to perform measurements at low frequencies with no bias to provide an additional observable 
of interest. 

Since electrostatic accelerometers suffer a bias, a technological upgrade has been proposed by 
Onera. It consists in adding to an electrostatic accelerometer a rotating platform which allows modu- 
lating the signal of interest and retrieving it without bias after post-processing. Using this principle, 
a measurement method and a post-processing method have been developed. The objective of this 
article is to validate these methods experimentally. To do so, a horizontally controlled pendulum was 
used to apply a known signal to an accelerometer mounted on a rotating platform. The processing 
of the experimental data demonstrates the ability to make acceleration measurements with no bias. 
In addition, the experimental precision on the unbiased acceleration obtained after post-processing 
corresponds to the precision predicted theoretically. 

Keywords Electrostatic accelerometer; Bias rejection; Modulation; Data processing; Precision; 
Colored noise. 

1 Introduction 

The Roadmap for Fundamental Physics in Space issued by the European Space Agency (ESA) in 2010 
put an emphasis on gravitation tests in the Solar System with missions to the outer planets. In this 
framework, it recommends the development of accelerometers compatible with spacecraft tracking at the 
10 pm.s~^ level at low frequencies. 

Indeed, interplanetary probes can be used as test masses whose trajectories are to be compared to 
theoretical predictions. Such a test has been performed by NASA with the Pioneer 10 and 11 probes: 
the outcome was a discrepancy with respect to the predictions of General Relativity O [3 HZ] . But these 
two probes lacked an instrument which could disentangle gravitational effects from non-gravitational 
ones. In order to improve the experiment made by the Pioneer probes, several space missions have been 
proposed jl] [SI HI [H [TH [23] , many of them embarking an accelerometer. The OSS mission [7^ relies, for its 
fundamental physics objectives, on the Gravity Advances Package jl6| . It is an instrument composed of 
an electrostatic accelerometer called MicroSTAR, based on Onera expertise in the field of accelerometry 
and gravimetry with CHAMP, GRACE, and GOCE missions [21], and a rotating platform, called Bias 
Rejection System (BRS). In orbit technology is used with technological upgrades aiming at reducing 
power consumption, size and mass. 

This instrument aims at measuring the non-gravitational acceleration of the spacecraft with a precision 
compatible with ESA requirement and with no bias. Thus, it provides an additional observable and allows 
removing, during the orbit restitution process, the effect of the non-gravitational forces on the trajectory. 
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enhancing deep space gravitation tests as well as gravity field recovery [13j. The measurements made by 
MicroSTAR or any other electrostatic accelerometer always include an intrinsic bias which may vary with 
time or temperature. It calls for an additional subsystem able to remove the bias through the modulation 
of the signal of interest. This is the role of the Bias Rejection System. In order to remove properly the 
bias of MicroSTAR, the Bias Rejection System rotates the accelerometer following a carefully designed 
periodical pattern [T3]. In terms of measurement noise, this operation selects the noise of MicroSTAR 
around the modulation frequency. After post-processing and for a modulation period of 10 min, it 
allows making absolute measurements with a white noise whose Power Spectrum Density (PSD) level is 
10^^*^ m.s^^.Hz^^/^ with a cut-off frequency equal to 8.3 x 10~^ Hz. This corresponds, for an integration 
time of 3 h, to a precision of 1 pm.s"^ and an exact measurement accuracy [T^. When taking into account 
the integration of the instrument in the spacecraft, and in particular the alignment accuracy (< 1 mrad), 
the positioning accuracy and the spacecraft self-gravity, a global precision of 10 pm.s~^ is expected. 

The spacecraft of the OSS mission has been specifically designed to reach this level of precision. The 
three main drivers taken into account are to (i) provide the lowest and most axisymmetrical gravitational 
field as viewed from GAP, (ii) make coincide the dry mass center of gravity, the propellant center of gravity, 
the radiation pressure force line and the GAP, and (iii) ensure a stable and reliable alignment between the 
GAP and the high gain antenna to ensure consistency between radio science data and non-gravitational 
acceleration measurements. This led to the architecture presented in [3: the platform is built as a fiat 
ring and the GAP is on a settable plate at the center, held by thermally stable struts. In addition, 
four Hydrazine tanks are distributed symmetrically around the GAP, with fiuidic interconnections. It 
allows balancing the quantity of propellant and thus controlling the position of the center of mass of the 
spacecraft. 

The goal of the research project presented in this article is to validate experimentally the data pro- 
cessing method developed in |15] . First, it will be shown that it is experimentally possible to completely 
remove the bias from the measurements. Then, the precision on the unbiased acceleration obtained after 
post-processing will be compared to the predicted precision given the noise of the devices used in this 
experiment. These results will demonstrate the possibility to make unbiased acceleration measurements 
in real conditions using in space technology. Therefore, they are a key milestone toward embarking the 
Gravity Advanced Package on an interplanetary mission. 

Since the Gravity Advanced Package is not available yet, commercial devices were used to mimic it: 
the "measurement device" which replaces the Gravity Advanced Package and the overall experimental 
setup is described in section [2] In section [3] the experimental methodology as well as the data processing 
is presented. Finally, the results are discussed in section |4] 



2 Experimental setup 

The overall experimental setup is shown in fig. [T] shows. The "measurement device", used to mimic 
the Gravity Advanced Package, is described in Section |2.1[ A pendulum controlled by an electrostatic 
accelerometer, called DM1, is used to impose a monitored external acceleration to the measurement device 



(cf. Section 2.2 1 



2.1 Measurement device 

The measurement device is composed of: 

• A rotating platform RGVIOOBL from Newport which acts as the Bias Rejection System. The 
axis of rotation is called x and is vertical. The rotation is parameterized by an angle called 9. The 
angle repeatability is equal to 0.0003°. Given the noise of the Q-Flex (cf. Fig. [3|, the rotating 
platform does not introduce uncertainty in the measurement process. 

• Two Q-Flex QA-700 accelerometers from Honeywell |10) . They are mono-axial accelerometers and 
are positioned orthogonally so as to have bi-axial measurements. The Q-Flex A and B measure 
local gravity along, respectively, the z and y axis of the frame attached to the rotating part of the 
rotating platform. 

As the output of the Q-Flex is a current proportional to the measured acceleration, the measurement 
chain of one Q-Flex accelerometer is composed of a current- voltage converter and a gain. In order to obtain 
the scale factor of the whole instrument, the calibration process is made in two steps. First, the scale 
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(a) The attitude of the pendulum is controlled by actua- 
tors (circled in red). The DM1 (circled in green) is used 
in the control loop to monitor the attitude of the pen- 
dulum. This set up allows removing the seismic noise. 
The two Q-Flex accelerometers (A and B) mounted on 
a Newport rotating platform (circled in blue) are used 
to perform the measurements in this experiment. The 
electronics are not visible on this picture. 




(b) The local reference frame {Xo,Yq,Zq) is defined such 
that Xo is vertical. The angles </> and a (cf. Section [2.2[ l 
measure the inclination of the pendulum with respect to 
the horizontal plane. 



Figure 1: Experimental setup 



Q-Flex + converter 0.8376 V.s^.m^^ 

gain 1010.9 

Total Ga = 846.73 V.s^.ni-i 



Table 1: Measured scale factor of the Q-Flex A 



factor of each accelerometer plus its current-voltage converter is computed. To do so, the measurement 
axis of the Q-Flex is set vertical so that the local gravity field g is measured ; then it is fiipped so as to 
measure —g. By subtracting and adding these two measurements and knowing the local gravity field, the 
gain and the bias of the Q-Flex with the current-voltage converter are known. Then, the scale factor of 
the gain is obtained by measuring the output (in V) knowing precisely the input (in V). The numerical 
values of the scale factors are summarized in tables [l] and |2| Ga and Gb are the experimental values 
whereas Ga and Gb are the exact ones, which are unknown in practice. Let introduce the measurement 
scale factors on each axis ky = {Gb — Gb)/Gb and kz = {Ga — Ga)/Ga, which are small compared to 1. 

The numerical values have been chosen such that the scale factors of the whole measurement chains 
for the Q-Flex are comparable to the ones of the DM1 (cf. next section). It is also possible to know 
approximately the numerical values of the bias of the whole measurement chains for the Q-Flex A and 
B, respectively 6a = -224.50 V and 6b = -52.19 V. 

Finally, data need to be digitized for post-processing. To do so, the low-pass filter SR640 |2T| is used 
to perform anti-aliasing before digitizing the data with the Analog-to-Digital Converter (ADC) NT 6033 
card lllj. The characterization of these devices in term of bias, scale factor and noise has been performed 
and is used in the data processing. For every experiment, the cut-off frequency is equal to the half of the 
sampling frequency in order to avoid aliasing [TH', p. 77]. The measurement chain is shown in fig. |2] 

In order to compare the predicted precision [15_ to the experimental one, it is necessary to characterize 

the noise of the measurement chain described above. Fig.|3]shows the measured Power Spectrum Density 

- 2 ~ 2 

(PSD). It can be approximated by the analytical functions S'q^(/) = S{f)/G^ and Sq^ = S{f)/G^ 

(in m^.s~"'.Hz~^) corresponding respectively to the Q-Flex A and B with 



S{f) = 3.1925 X 10-^ V^.Hz-i + " ^'^^ (1) 



9.052 X 10-« V^Hz' 



2 TT^O.6384 



1.6384 
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Q-Flex - 


h converter 0.8120 V.s^.ni^^ 


gain 


990.54 


Total 


Gb = 804.32 V.s^.m-i 



Table 2: Measured scale factor of the Q-Flex B 




- Data recording 

- Data processing 

- Control of the BRS 



Figure 2: Schematic of the measurement chain. The computer is used for data recording and processing 
and for controlling the rotating platform. A low-pass filter is used before the Analog-to-Digital Converter 
(ADC) to avoid aliasing. 



For frequencies above 2 x 10^^ Hz, the analytical model does not fit well with the measured PSD. 
Indeed, at these frequencies, it is not the electronic noise of the Q-Flex which is predominant but the 
seismic noise filtered by the pendulum response function. It is however not a concern because these 
frequencies are too high to play a role in this experiment. For frequencies below 2 x 10^^ Hz, the noise of 
the Q-Flex is far larger than the noise of the pendulum (cf. Fig.|4]). This means that in this experiment, 
only the noise of the Q-Flex accelerometers plays a role. 

2.2 Horizontally controlled Pendulum 

The pendulum is used as a mean to deliver a known signal to the measurement device with a noise 
characterized by fig. |4j To achieve this level of noise, it is mounted on concrete decoupled from the 
ground. Moreover, the outputs of an electrostatic accelerometer, called DM1, are used as sensors: the 
servo-loop is designed such that the outputs in the horizontal plane {Y and Z axis of the DM1) are equal 
to zero. This accelerometer is an engineering model of the Aristoteles mission [20_. Tabled gives the 
experimental values of the scale factors of the two horizontal axis of the DM1, Gy and Gz- The exact 
values are named Gy and Gz- As for the Q-Flex, let's also introduce the measurement scale factors on 
each axis ky — [Gy — Gy)/Gy and kz — [Gz — Gz)IGz- 

Y axis Gy = 805 V.s^.m-^ 
Z axis Gz = 825 V.s^.m-^ 



Table 3: Scale factors of the two horizontal axis of the DM1 accelerometer [5]. 

Because of the bias hy and hz of the DM1 accelerometer along the Y and Z axis, the horizontal axis 
of the DM1 are perpendicular to local gravity with an offset. In addition, it is possible to incline the 
DM1 by a known angle. To do so, secondary entries on the Y and Z axis of the DM1, called Vy and Vz 
(in V), allow to put offsets. Given the scale factors of Table |3] the inclination angle and a of the DM1, 
respectively around the axis Y and Z, are equal to 



a 



bz-Vz 
gGz 

by - Vy 

gGy 



1 

g 
1 

9 



bz^ 
Gz 

bY_ 
Gy 



1 

l + kz 

1 

1 + fcv 



cz 



-Cy 



(2a) 
(2b) 
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Experiment 

- - - Model 




10"' 10"' 10"' 

Frequency (Hz) 

Figure 3: Square-root of the Power Spectrum Density (PSD) of the noise of the Q-Flex A and its 
measurement chain. The analytical model is given in equation (nl. 




Figure 4: Square-root of the Power Spectrum Density of the pendulum noise. Below 2 X 10^1 Hz. the 
pendulum is servo-looped on the outputs of the DM1 accelerometer and the noise level goes down to 
approximately 3 x 10~* m.s""^.Hz""^/^. Above 2 x 10~^ Hz, the pendulum is servo-looped on the ground 
and the seismic noise, filtered by the pendulum response function, is dominant. 



where cy = Vy/Gy and cz — Vz/Gz are offsets in term of acceleration. A typical amplitude for Vy and 
Vz is 1.5 V, which leads to an inclination of the order of 2 x 10^"' rad and a projection of local gravity 
in the Y-Z plane of the order of 2 x 10""^ m.s~^. 



3 Methods 

3.1 Modeling of the experiment 

As explained, the attitude of the DM1 with respect to the local reference frame is parameterized by the 
angles a and cj). Because of the experimental setup, the alignment between the DM1 and the rotating 
platform is not perfectly known. Therefore, the attitude of the rotating platform with respect to the 
local reference frame is parameterized by the angles d = a + aa and = + 0Oi where ao and 0o £^re 
constant but unknown. For these experiments, a has been set to 0, which means that only rotations of 
the pendulum around the axis Y have been applied. 

Calling rriz and niy the measurements (in m.s^^) made respectively by the Q-Flex A and B, they are 
equal at first order to 

f my — ky [ay cos{6) + az sin(6')] + by (3a) 
\ rriz = kz [-ay sin{9) + az cos{9)\ + (3b) 
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Figure 5: Outputs of the Q-Flex as a function of tlie angle 9 ioi a — (f) — 0. The amplitude of the curves 
is linked to the external acceleration and the offset of the curves is due to the bias and the wobble of the 
measurement device. 



with ay and az the external accelerations (in m.s ^) 

ay ~ gao (4a) 
az = -gip (4b) 

and by and the measurement bias (in m.s~^) of the whole "measurement device" (Q-Flex + rotating 
platform) : 

f by = bs/GB + gPs (5a) 
1 - bA/GA - gfiA (5b) 

where Pa and Pb are the wobbles of the Q-Flex A and B on the rotating platform. The wobble is defined 
as the angle between the measurement axis of the Q-Flex and the plane perpendicular to the rotation 
axis of the rotating platform. The bias 6^ and bs are known experimentally (cf. Section 2.11, but 
and Pb are not, such that the values of by and bz are not known a priori. 

Equations ([3| are identical to those studied in [15], which considered the measurement made by 
the Gravity Advanced Package (GAP) on an interplanetary spacecraft. Therefore, this experiment is 
representative of the measurements made by the GAP during a space mission and the conclusion drawn 
in this article can be applied to the space experiment, except that the noise of the Q-Flex is higher than 
the noise of the GAP. Fig|5]shows the variation of my and with the angle controlled by the rotating 
platform. The data are fitted, and it leads to the following numerical values: 

my = [-2.94cos(6i) - 2.98sin(6l) - 14.33] x 10~* m.s"^ (6a) 
rriz = [-3.34 cos(6l) + 3.02 sin(6l) + 10.54] x 10"'' m.s^^ (g^^) 

Using the offset of the two curves and knowing the bias of each Q-Flex, it is possible to determine the 
values of the wobbles : Pa = -8.07 x 10"^ rad and Pb = 1-40 x IQ-^ rad. 

3.2 Description of the measurements 

Concerning the applied external acceleration on the Z axis, cz, controlled by the inclination cj) of the 
pendulum, two types of signals have been used. In one case, a constant inclination of the pendulum has 
been applied with different values for the offset Vz of the DM1 accelerometer: V, 0.5 V, 1 V and 1.5 \Q 



The data gathered with these signals are used in Section 4.3 In the second case, the inclination of the 
pendulum had a sinusoidal and a triangular variation with an amplitude of 1 V and different periods: 
0.75 min. 1 min, 5 min and 60 min. These data are use in Section 14.21 



respectively. 
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Figure 6: Calibration signals 9{t) used for the experiments. Two periods are represented, separated by 
circles (o). The total time spent rotating is always equal to 20 seconds per period. The signal 6(a) and 
|6(b)| are exactly similar except that the origin of the period is different, but they do not have the same 
properties as far as post-processing is concerned (cf. Section 4.2 1. The signal 6(c) is different from the 
signal [6(a) because it does not display a periodicity of 0.5 arbitrary unit. 



The time patterns used for the angle 9, which parameterizes the rotation of the Q-Flex, are the one 
introduced in [15 . They are represented in fig. [6j Several modulation periods have been used: 1 min, 
5 min, 20 min and 50 min. The characteristics of each signal in term of ability to remove the bias from 
the data is discussed in Section 14.21 

The sampling frequency used was mainly 10 Hz. Other sampling frequencies have been used so as to 
demonstrate that the precision on the post-processed quantities do not vary with this parameter under 
some conditions. 

For all the data used in this article, the duration of the measurement was equal to 100 times the 
modulation period (e.g., for a modulation period of 50 min, the duration of the experiment was 3 days 
and 11 hours). 



3.3 Data processing 

Ac = diag[cos(6lfc)] and A, = diag[sin(6'fc)], ^ € ||1; 7V||. (7) 

When the rotating platform is moving, it may induce additional signal due to off-centering as well as 
vibration. Therefore, only the measurement made when is constant are considered. Given the signals 
offig.|6] A, =0. 

Because the scale factors ky and kz are unknowns, it is not possible to determine the quantities aY 
and az- Therefore, the quantities considered in the rest of this articles will be ay = (1 + ky)aY and 
az = {l + kz)az- Considering these quantities is natural because the values for the angle 9 are and 180 
degrees (cf. Fig. |6|. 

The overall goal of data processing is to find the projection of the vectors eIy and az on a vector 
subspace (of dimension pa < N) whose basis is made of the column of a matrix Va € M N^p^ (M) , which are 
supposed to be orthogonal for the usual scalar product on M^. As a result, the goal is to find the numerical 
values of {VaVa)~^Va' Sly and {VaVa)^^Va'az knowing niy and niz {M' is the matrix transpose of M). 



In this article, the choice of Va (cf. Section 4.2 1 will allow retrieving the mean value of the acceleration 
without bias and/or the slope of the acceleration over one modulation period. 

It has been demonstrated [TS] that it is possible to remove perfectly the bias from the measurements 
and that 

K'av = Va'Acmy (8a) 
K'az = Va'Acin^ (8b) 
under the following conditions 

K'A,b, = 0, with K e {y; z}. (9) 

In the same way, it is possible to retrieve the bias of the instrument. As for the acceleration, the goal 
is to find the projection of the vectors by and on a vector subspace (of dimension pb < N) whose 
basis is made of the column of a matrix Vf, S A^Ar_p^(M), which are suppose to be orthogonal for the usual 
scalar product on M^. Under the following conditions 



Vb'Aca„ = 0, with k G {Y; Z}. 



(10) 
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it is possible to retrieve the bias: 

Vb'hy = Vb'niy (11a) 
Vb'K = Vb'my (lib) 

As for the external acceleration, VJ, will be chosen such that the mean value of the bias and/or the slope 
of the bias over one modulation period are retrieved after post-processing. 

It is a priori not possible to know if conditions Q and ( 10 1 are fulfilled since the temporal evolution 



of the bias and of the external acceleration may not be controlled. In this experiment however, it is 
possible to assume that the bias belongs to the subspace generated by the columns of a matrix Vt and 
that the external acceleration belongs to the subspace generated by the columns of a matrix Va- Under 



these assumptions and given that = 0, the previous conditions (|9| and ( 10 1 come down to Va'AcVb = 
and Vb'AcVa = 0. 

Finally, it is essential to notice that the post-processed acceleration and bias are obtained separately 
and under different conditions. For example, it is possible to retrieve correctly the external acceleration 
without bias if Vj AcVb = whereas the post-processed bias will be mixed with the external signal because 

Vb'AcVa ^ 0. 



4 Results and discussion 



Three goals are pursued with the experimental setup described above. First, a validation of the modeling 
made in Section |3.1| is presented. Second, it will be shown in Section |4.2| that obtaining unbiased mea- 
surement is possible applying the post-processing scheme described above. Finally, in Section |4.3[ the 
experimental uncertainty of the unbiased measurements will be shown to follow the theoretical prediction. 



4.1 Calibration of the relative orientation 

It is interesting to aggregate all the experimental data in order to calibrate the unknown scale factors. 
For each condition with a constant inclination of the pendulum cz, an experimental value for az has been 



computed from the data by applying equation ( 8b I . Considering equations ( 2a I and ( 4b I , the theoretical 
link between these two quantities is 



1 



-cz - (1 



bz_ 
Gz 



1 + kz 

Experimentally, the following numerical values can be computed from the data 

= 0.964 XCZ + 2.8757 x lO"'^ m.s"^. 



(12) 



(13) 



The experimental value of (1 -I- fcz)/(l + kz), which is close to 1, demonstrates that the numerical values 
Ga and Gz used to process the data are close to the real ones. Unfortunately, this fit does not allow to 
compute independently the value of bz and (j)Q. 



4.2 Experimental validation of the demodulation process 



In Section 3.3 the equations i\9[ and (10 1 gave the conditions under which it is possible to retrieve the 



unbiased acceleration and the bias of the instrument from the data. To be more specific, let's introduce 
the following matrices: 



Vi = 







(14) 



and 



V2 



(15) 



where Ig is a matrix of A4q,i{^) whose coefficients are 1, and tg is a matrix of Aiq^i(U.) such that 
tq^ = (fc — q/2)6t. q is the number of data points during one period of the calibration signal. Assuming 
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that a vector x belongs to the subspace generated by the matrix V2 means that it is an afEne function of 
time on each modulation period. For a signal whose characteristic duration of variation is long compared 
to the period of the calibration signal, it is legitimate to assume that this signal can be approximated 
locally by an affine function and therefore that it belongs to the subspace generated by the matrix V2 (it 
will be done below with a sinusoid). 



For the calibration signal of fig. 6(a) ViAc^l — and ViAcV2 ^ 0. It means, for example, that it is 
not possible to recover correctly the mean bias over a modulation period if the external acceleration 

(D- 



is an affine function of time on each modulation period (cf. Fig. 



For the calibration signal of fig. |6(b)[ ViAcV2 = and V2ACV2 ^ 0. Fig. |8(b)] provides an illustration 
of the case. 



For the calibration signal of fig. 6(c) V2ACV2 = 0. In this case, it is possible to recover correctly 
the mean and the slope of the unbiased acceleration and of the bias over a modulation period under 
the assumption that they are affine functions of time on each period of the calibration signal. 



4.2.1 Perfect demodulation 

First, let consider a case for which V2ACV2 = 0. To do so, the external signal is a triangular variation and 
the modulation signal is the one given in fig. |6(c) The results of the data processing are plotted in fig. [7j 
Several features shows that post-processing allows recovering correctly the mean and the slope of the 
acceleration and bias over one period of the calibration signal. First the post-processed mean bias is 
constant, as expected since the temperature was constant (cf. fig. 7(a)[ ). Second, the post-processed 
mean acceleration precisely corresponds to the external one with no bias, as shown by fig. |7(b)[ Finally, 
fig. 7(c) shows that it is also possible to retrieve the slope of the acceleration. In this figure however, it is 
possible to see that the slope of the bias has a variation which looks like the one of the external signal: it 



means that the conditions (10 1 are not perfectly verified. This may come from the fact that according to 
fig. 7(c)| the slope of the external acceleration is not perfectly constant but has random variations which 
interfere. 



4.2.2 Interference of the signal and the bias during post-processing 



To illustrate the behavior of post-processing when the conditions ([9| and (10 1 are not respected, let's 
consider the following experimental conditions: the external signal is a sinusoidal variation with a period 



of 60 min with a peak-to-peak value for Vz of 1 V; the calibration signal is the signal of fig. 6(a) with a 
period of 5 min. 

The data are processed in two different ways. In one case, the signal of fig. 6(a) is used, and in the 



other case, the period frame is shifted by 1/4 of a period and the data are processed using the signal of 
fig. |6(b)[ The results are plotted in fig. [s] 

Because the sinusoid cannot be considered as a constant function on each modulation period, Vi A^a^ ^ 



for the signal 6(a) As a result, the post-processed mean bias is the addition of the real mean bias and 
an additional term, which is {Vi'Vi)^^Vi' AcSlz- This is a column vector and each value correspond to a 
sampling point. It can be approximated by the value of the following function computed at the sampling 
points: 



h{t) 



t+T/2 



c% cos{ujX)dX 



t+T 



t + T/2 



c% cos(wA)dA 



24 



cos 



\2 



sin iot 



LUT\ 



(16) 



where r is the modulation period, uj is the pulsation of the sinusoid and c% its amplitude. For the 



experimental conditions of fig. 8(a) the peak-to-peak value of h is 3.15 x 10 ^ m.s ^ and its time shift 



is 17.5 min. It corresponds perfectly to the sinusoidal part of the post-processed bias of fig. |8(a) 



The result of the data processing with the calibration signal [6(b)] are plotted in fig. |8(b)[ In this case, 
because ViAca^ = for the the signal [6(b)] there is no mix of the external acceleration with the bias: it is 
possible to retrieve the correct values of the mean bias and of the mean acceleration after post-processing. 
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(a) Post-processed measurements giving tlie mean ac- 
celeration and the mean bias over a period of the cal- 
ibration signal. Equation \12\ was used to correct for 
the unknown alignment coefficients. 
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(b) Difference between the post-processed mean accel- 
eration and the real mean acceleration. 
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(c) Post-processed measurements giving the slope of the 
acceleration and of the bias over a period of the calibra- 
tion signal 



Figure 7: Experimental conditions: the external signal is a triangular variation with a period of 60 min 
and an amplitude Vz — 1 V] the calibration signal is the signal of fig. 6(c) with a period of 1 min. In 
fig. (a) and 1(c) I the line (-) represent the known external signal, the dots (•) represents the post-processed 
acceleration, and the circles (o) represented the post-processed bias. 



4.3 Uncertainty on the measurements: comparison to theoretical prediction 

One of the important goal of this experiment is to be able to verify experimentally the precision on the 
mean unbiased acceleration. To do so, data acquired with a constant inclination of the pendulum are 
used and aggregated. The theoretical precision on the mean unbiased acceleration for an integration time 
T measured after data processing can be expressed using the PSD of the Q-Flex A noise (cf. eq. (jlj) and 



the calibration signal. In this section, the focus will be on the calibration signal 6(c) because it is the 
most appropriate for our purpose, as demonstrated in the previous section. The complete formula giving 
the uncertainty can be found in [15_ and is used for data processing. For simplicity however, let us recall 
the simplified expression giving the uncertainty a for an integration time T and a modulation period t: 

a{T,T)^JlsQ(^\ (17) 



The general formula is an integral on the frequency of the noise PSD multiplied by the square of the norm 
of the Discrete Time Fourier Transform (DTFT) of the calibration signal. Moreover, since it has been 
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(a) Post-processed mean quantities over a period of the 
calibration signal. The data processing was made using 
the calibration signal of Fig. |6(a)| The peak-to-peak 
value of the variation of the post-processed bias is equal 
to 3.06 X 10-* m.s-2. 
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(b) Post-processed mean quantities over a period of the 
calibration signal. The data processing was made us- 
ing the calibration signal of Fig. |6(b)| The bias is not 
perfectly constant because the approximation that the 
sinusoid is afRne with time on each modulation period 
is not perfectly correct, especially at the extrema. 



Figure 8: Experimental conditions: external signal is a sinusoidal variation with a period of 60 min with 
an amplitude = 1 V; the calibration signal is the signal of Fig. 6(a) with a period of 5 min. The 
line (-) represents the known external signal, the dots (•) represents the post-processed acceleration, and 
the circles (o) represented the post-processed bias. Equation ( 12 1 was used to correct for the unknown 
alignment coefficients. 



shown [15 that the noise on post-processed acceleration is a white noise, (t(t, Ti) = y/T^jTi x ^(r, T2). 
Therefore, the only integration time considered in the following is 1 hour. 

This theoretical value for the uncertainty is compared to the experimental data in fig. [9j The exper- 
imental data are obtained by computing the standard deviation of the experimental values of az when 
the pendulum inclination is constant. Since the noise of the pendulum is much lower than the one of the 
Q-Flex, the measured values for az changes from one modulation period to the other because of the mea- 
surement noise of the Q-Flex. By computing the standard deviation, one has access to the experimental 
precision on az- 

Fig. [9] shows that the experimental values are in very good agreement with the theoretical prediction 
over a range of modulation period which is considered for the Gravity Advanced Package. Therefore, 
these experiments validate the theoretical approach developed in a previous article [15^ as well as the 



theoretical formula (17). The fact that the Q-Flex have a noise level much more higher than the noise 



level of MicroSTAR does not jeopardize this experimental validation of principle. 



Frequency (Hz) 


Precision (m.s ^) 


10 


1.6631 X IQ-'^ 


100 


1.5718 X lO"'^ 


200 


1.2524 X 10-^ 



Table 4: Influence of the sampling frequency over the precision of the unbiased acceleration. The precisions 
are given for an integration time of 1 hour. Experimental conditions: external signal is a constant with 
= V; the calibration signal is the signal 6(c) with a period of 1 min. The theoretical precision is 
1.8073 X 10~^ m.s^^ 



The effect of the sampling frequency has also been investigated. Because the harmonics of the calibra- 
tion signal in the frequency domain decreases rapidly, it is expected that the precision does not depend on 
the sampling frequency as far as it is two order of magnitude larger than the modulation frequency. That 



is why the sampling frequency does not appear in the approximated formula (17). In these experiments, 
the larger modulation frequency is 1/60 = 0.0167 Hz and the explored sampling frequencies are 10, 100 
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Figure 9: Precision (Icr) of the mean unbiased acceleration as a function of the period of the cahbration 
signal for an integration time of one hour. The bars show the confidence interval at 99 %. The data have 
been obtained with a sampling frequency of 10 Hz and with the calibration signal described in fig. |6(c) 
Each point comes from the statistic over 200 to 700 data points depending on the modulation period. 
The theoretical uncertainty is also plotted using the general formula given in |15j . 



and 200 Hz, which are two orders of magnitude larger than 0.0167 Hz. The results are summarized in 
Table [4] and show that the precision does not depend significantly on the sampling frequency for the range 
explored in these experiments. 



5 Conclusion 

The set of experiments presented in this article was designed to validate the data processing scheme 
developed in order to remove the bias from the measurements made with an electrostatic accelerometer. 
To do so, two Q-Flex accelerometers mounted on a rotating platform were placed on a pendulum whose 
inclination was precisely controlled via a high-precision electrostatic accelerometer. 

The data gathered provided an experimental validation of the method. First, the ability to recover, 
after post-processing, the mean and the slope of the bias and of the unbiased acceleration has been 
demonstrated with experimental data. Then, concerning the precision of the unbiased acceleration mea- 
surement, the experimental precisions have been compared to the theoretical predictions and a solid 
agreement has been obtained. 

The experimental validation of this method devised to make unbiased acceleration measurement with 
an electrostatic accelerometer opens new opportunities. It will allow improved orbit restitution of inter- 
planetary probes using an instrument with a strong heritage. 
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